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Abstract 

We investigate the wicking in granular media by considering layers of 
grains at the surface of a liquid, and discuss the critical contact angle below 
which spontaneous impregnation takes place. This angle is found to be 
on the order of 55° for monodisperse layers, significantly smaller than 90°, 
the threshold value for penetrating assemblies of tubes: owing to geometry, 
impregnating grains is more "demanding" than impregnating tubes. We also 
consider the additional effects of polydispersity and pressure on this wetting 
transition and discuss the corresponding shift observed on the critical contact 
angle. 

Introduction 

spontaneous invasion of pores by liquids is driven by a reduction of surface energy. 
This phenomenon is relevant to many natural applications such as soil imbibition 
by rain [[ll|2l[3||, soil decontamination [4J or plant physiology [|5l[6]|. It is also 
of broad industrial relevance for oil extraction, civil engineering [[Z||, absorbent 
consumer goods [|8^ i9J, paper and textile design [|10l [U, chromatographic [13 
and microfluidic processes [13 J. An important application in the food industry 
is the reconstitution of a beverage upon wetting of a dehydrated powder |[T4| . 
Depending on the field of interest, impregnation can either be beneficial (food 
industry) or detrimental (civil engineering). 

Figure 1 compares wicking in a porous media made of spherical beads and in a 
capillary tube. For large contact angles, there is no wicking (Figure la). Wicking is 
observed when the contact angle is below a critical value: ninety degrees for the 
capillary tube (Figure lb) and a much lower one for the beads (Figure Ic). This 
experiment emphasizes the role of the geometry in wicking: a porous medium 
made of grains or beads differs substantially from cylindrical capillaries, because 
cross sectional area and wall orientation vary along the pores. This modifies the 
condition of wicking, as also shown for other special porous media [[T5l[T6l[T7l| . 
However, the classical theories for cylindrical capillary tubes |T8ilT9ll20B are still 
broadly applied to powders. 
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Figure 1: A silanised capillary tube and a porous media made of silanised glass 
beads are put in contact with baths of various water-ethanol mixtures, changing 
the contact angle 6. (a) No wicking is observed at high contact angles, (b) Wicking 
in the tube occurs when the contact angle is below 90^ but the grains remain dry. (c) 
Spontaneous impregnation in the porous media is observed only below a critical 
contact angle 6^ significantly smaller than 90^. Scale bars correspond to 5 mm. 



In the general case of a compact pile of grains, the maximum penetration 
height can be computed by adapting appropriately Jurin's theory [16J. In the 
same spirit, the classical kinetics of Lucas-Washburn [|2lH22|| can be modified by 
considering the effect of the variable cross section of the granular pores on the 
viscous dissipation, which generates the rich impregnation kinetics of layered 
beads [^^. In these works, the capillary driving force is often considered as 
constant along pore length. Reyssat et al. [24J considers for instance complete 
wetting and includes an adjustable parameter in the capillary term. Fries assumes 
that spontaneous impregnation in granular beds occurs if the contact angle of the 
liquid with the channels is lower than ninety degrees, that is the condition for 
pores having walls parallel to the pore axis. Tsori [25J and Czachor 112611 studied 
theoretically the capillary penetration in channels with varying cross section, 
including sinusoidal capillaries. The former showed the existence of multiple 
equilibrium positions, while the latter predicted the existence of a critical wetting 
angle for capillary rise that depends on the sinusoidal wall waviness. Lago et al. 
[|27| investigated the optimal wetting angle for the highest capillary rise, that is 
15.4"^ due to the geometry of a compact pile of spherical grains. Ban et al. [I28| 
and Shirtcliffe et al. [29J discussed spontaneous wicking of liquid in such piles. 
They predicted the value of a critical contact angle 6^ < 90^ for the wicking in this 
geometry, and produced experimental evidences of this result. These works will 
be further commented in section 3 and in section 4. 

In this article, we report experiments and discuss a model, showing the exis- 
tence of a critical equilibrium contact angle above which the impregnation of a 
heap of grains can be blocked, as seen in figure 1. We also discuss how this critical 
contact angle depends on polydispersity and hydrostatic forcing. These results 
can be used either to improve or to block capillary invasion, depending on the 
system of interest. 
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1 Experimental set-up and protocol 




Figure 2: (a) Experimental set-up and illustration of the impregnation criterion, (b) 
Glass beads with a radius Ri = 52 ± 2 ^m, as observed at the surface of a liquid 
bath, (c) Corresponding distribution of beads radii for a sample of 100 beads. 

In order to study the impregnation of a granular media, we deposit with a pa- 
per cone a dry pile of grains at the surface of a liquid bath filling a centimetric 
transparent cell (Fig. |2^). The liquid-grain interface is observed with a camera, 
magnification x2. Our wicking criterion is the detachment of grains from the 
interface, which then fall in the bath due to a density ps = 2450 kg/m^ higher than 
the one of water. Conversely, the liquid does not penetrate the pile when all the 
grains remain at the surface. 

The grains are borosilicate glass beads (from Sigmund Lindner), sieved in 
order to reduce their polydispersity (Figures |2j:) and [2]:). Several radii R are used: 
25 ± 2 jim, 52 ± 2 ^m, 100 ± 5 }im and 256 ± 13 }im. The glass beads, initially 
hydrophilic, are silanized with 1H-1H-2H-2H perfluoro-octyltriethoxysilane to be 
hydrophobic, following the protocol by Qian et al. [30J. 

The wetting properties of the system are tuned by using mixtures of water and 
ethanol. The ethanol volume fraction i;^^ in water is controlled by measuring the 
density of the mixture, referring to tables in Handbook of Chemistry and Physics 

[|3T1l . The capillary length a = (where 7 and p are the surface tension and 

density of the liquid, and g the gravity constant) varies from 2.8 mm for deionized 
water (i;^^ = 0) to 1.7 mm for pure ethanol {vf = 1), always larger than the beads 
radius. 

In order to measure the contact angle 6, single beads are placed on the surface 
of the liquid (Fig. [3]), and 6 is deduced from the distance S = R cos 6 between 
the bead center and the liquid interface (for R <^ a). Pictures are taken with a 
Ricoh GX200 camera trough a binocular x20 (Fig. [3]). The position of the liquid-gas 
interface is determined by the reflection of the beads on the surface. Measurement 
of contact angle for each i;^^ is repeated on 10 different beads. The evolution is 
presented in Fig. [4| The contact angle 6 with deionized water is 105^ ± 5^, and it 
falls down to 35"^ ± 5"^ for pure ethanol. Mixtures of water and ethanol allow us to 
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(a) 17^ = 0,^ = 105° 



(b) Vf = 0.55,0 = 61' 



(c) Vf = 0.95,0 = 35"" 



Figure 3: Measurement of the contact angle between a glass bead (here, R = 
256 fim) and a mixture of water and ethanol of ethanol volume fraction Vf. The 
dashed line indicates the liquid/ air interface. 

obtain intermediate contact angles, which continuously decrease from 105"^ to 35^ 
asvf passes from to 1. 

2 Experimental results 






(a) Vf = 0.4,0 = 75'' ±5"" 



(b) 17^ = 0.65, = 58° ±7° 



(c) 17^ = 0.7,0 = 56° ±5° 



Figure 5: Snapshots of impregnation experiments for different ethanol volume 
fractions Vf. A pile of grains is made at the interface between the bath and air. 
The apparent thickness (in black) of the interface comes from the menisci on the 
edges. Wicking is deduced from the presence of grains falling in the bath, (a) no 
impregnation, (b) limit case, (c) impregnation. Glass beads radius is 256 ± 13 }im 
for each experiment. 



We deposit a uniform layer of dry grains at the surface of a liquid bath of ethanol 
fraction Vf (Fig. [5]). A meniscus forms along the cell walls, which generates the 
dark zone observed in Fig. [5} The actual thickness of the dry powder represents 
only a few layers of beads (typically 4). Two regimes are observed. In Fig. |5^, 
no bead falls into the bath and the pile remains dry. In Fig. [5]:, beads detach 
and fall as the liquid invades the pile. Fig. |5]3 shows the limit of impregnation. 
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Figure 4: Contact angle 6 between glass beads of radius R = 52 }im (red circles) 
and R = 256 }im (blue triangles), and water-ethanol mixtures, as a function of the 
ethanol volume fraction Vf. The contact angle can be varied continuously between 
105^ ± 5^ and 35° ± 5° by increasing the ethanol fraction. 



where only very few beads fall. The wicking transition occurs for an ethanol 
fraction i;^ = 0.65 ± 0.03 for R = 256 ^m, which corresponds to a contact angle 
6^ = 58° ±7°. ^'^^ is a critical angle: if we have Q < 6^, wicking is observed; if not, 
the pile remains dry. The results for different grain radii are presented in table [ij 
showing that 9^ hardly depends on R when grains are smaller than the capillary 
length {R << a). The highest value of the critical contact angle corresponds to 
the largest radius. Since the density ps of the grain is significantly higher than the 
density of the liquid, gravitational corrections are expected for the experiments 
with R = 256 }im. 



3 Model 

3.1 Interface equilibrium 

In order to explain these experiments, we first consider a single grain at the 
interface: since the spheres are denser than the liquid, gravity is balanced by 
interfacial forces. This is no longer possible if the sphere is surrounded by liquid, 
in which case it sinks into the liquid. 

Let 7 be the surface tension of the liquid, xp the angle between the equatorial 
plane and the radius that connects to the contact line, Rc the curvature of the 
interface and ^Rc the length of the liquid-air interface that connects the bead to 
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R {}im) 




r (°) 


15 ±1 


0.73 ± 0.01 


53 ±6 


51 ±1 


0.73 ± 0.01 


53 ±6 


100 ±5 


0.70 ± 0.01 


55 ±6 


256 ± 13 


0.65 ± 0.03 


58 ±7 



Table 1: Experimental results of the wetting transition in terms of critical ethanol 
fraction vJ, and corresponding contact angle 6^ , for different beads with radius R 
smaller than the capillary length a. 



the flat bath (Fig. [6]). The contact angle is satisfied at lower scales, thus the angle 
between the sphere and the interface must be 6. All these angles are linked by the 
geometrical relationship: 

^ + ^-iS=f. (1) 

Keller Il32ll has shown that the vertical projection of the pressure forces is equal 
to the weight of the volume of liquid bounded by the horizontal free surface of 
the bath, the wetted surface of the body, and the vertical cylinder of radius R cos \p 
and height ^ = Rc{l — cos j6) (Fig. [6]). At equilibrium, the vertical projection of 
forces on the body can be written as a balance between the surface tension force 
27zRj cos ip sin ^ and an effective weight ^nR^gp^ff, gathering the buoyancy and 
the weight of the sphere. This leads to: 

sin«=?^?^ = M^!^, (2) 

37 COS t/^ p a^Scosxp' 



where a=./J-is the capillary length. In our case, the radii are much smaller than 




Figure 6: Sphere at a curved interface. The vertical component of the buoyancy 
corresponds to the weight of fluid of the volume (in white, below the dashed 
line), that is the volume of the solid below the contact line, plus the one of the 
cylinder of radius R cos xp and height ^. It is the opposite of the weight of all the 
water displaced by the presence of the body minus the contribution of the surface 
tension (in gray). 
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the capillary length, and the ratio of densities is close to 1. Thus sin j6 is small, on 
the order of lO'^ for R ~ 100 }im. For 6 > 0, we can thus neglect the effect of the 
weight and consider a flat interface as indeed observed in the experiments (Fig. |3]). 
This corresponds toip + 6 = The sphere stabilizes at a height 5 below the flat 
interface (Fig. |7^), where S is given by: 

S = Rcose. (3) 




Figure 7: (a) Equilibrium state of a small bead at a liquid-gas interface, (b) Limit 
case for the impregnation in a 2-D close packing: the interface is tangent to the 
upper spheres. 



3.2 2-D pile 

Based on the experimental observation presented in Fig. [3]3), we now consider 
a compact pile of grains (Fig. [^). As the equilibrium discussed previously is 
achieved for the first layer, the liquid may reach the second layer of grains. If it 
happens, the contact line is no longer at equilibrium and it will move up into the 
pile until it satisfies both a flat interface and a contact angle 6 on the second layer. 
As it did for the first layer, it reaches the next spheres, and the same mechanism 
applies: no equilibrium is possible as long as there are dry grains. The liquid 
impregnates the whole pile, until only one monolayer of grains remains at the 
interface. As the liquid rises, the lower spheres get surrounded by liquid and thus 
detach from the rest of the pile and fall into the bath, as seen experimentally (Fig. 
|5]:). Conversely, if the liquid does not reach the second layer of grains, only the 
first layer is in contact with the liquid and no grain detaches from the interface. 

The limit between these two situations corresponds to 6 = ^20^ which 
the interface is tangent to the spheres of the second layer (Fig. [^). The distance 
between two equatorial planes is equal to R + S^ and it can also be written 2R cos a:, 
where oc is the angle between the line linking the sphere centers from different 
layers with the normal to the first layer (Fig. [^), thus giving: 

R + S"" = 2Rcos oc. (4) 
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In the compact 2-D case, the centers of monodisperse grains form equilateral 
triangles, which yields oc = 30°. This result, associated with equations ^ and Q, 
provides the value of the critical contact angle ^Jd impregnating a 2-D pile 
(assembly of infinite cylinders): 

= arccos(\/3 - 1) ^ 43^ (5) 

If the contact angle 6 is higher than this critical value (and smaller than 90"^), 
there is a local minimum of the surface energy of the system: impregnation can 
be blocked. The system may not reach its global minimum of energy, owing the 
existence of an energy barrier AE (per unit length) given by: 

AE sin e - sin 0^ + (0^ - 6) cos 6 

~ / (6) 



2nRj TC 
This expression becomes, at leading order inO — 



2nRy In 

For a bead's diameter IR ^ lOOjitn, the order of magnitude of the energy barrier is 
typically 2RAE ^ 10 J, much larger of course than thermal energy. 



3.3 3-D geometry 




Figure 8: Projection along the plane given by a median of the base of the tetrahe- 
dron and its summit for a 3-D pile. The dashed circle represents the two other 
spheres forming the base of the tetrahedron, out of this plane. 

Similar effects are expected in 3-D, and Shirtcliffe and al. [|29ll used surface 
free energies to predict the existence of a critical contact angle of wicking whose 
value is 50.73^ for a compact pile. We can obtain this result using geometrical 
considerations. In a 3-dimensional compact pile of grains, the centers of the 
spheres form a tetrahedral network. The discussion above is still correct as well as 
equations ^ and Q, but the relative position of the successive layers is slightly 
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modified, leading to a different value of cc, now given by sin a: = ^ (Fig. 
provides the critical angle for impregnating a 3-D pile: 

= arccos(A/ ^ -1)^ 51^ (8) 
V ^ 

Ban et al. [28] and Shirtcliffe et al. [|29ll derived the same critical condition for 
impregnation. Our experiments here are conducted differently: we measure the 
contact angle directly on the grains, as discussed in section 1. Moreover, we bring 
the grains in contact with the liquid without confinement. Using this experimental 
protocol allows us to detect impregnation as it occurs over a single grain layer, 
improving the precision of the measurement. 

The observed critical contact angles (Table [l]) are indeed close to 6^ ^ 51^, yet 
slightly larger. An additional deviation appears when R is higher than 100 /^m, re- 
vealing the influence of gravity, neglected in the model above. Different hypothesis 
can be proposed to explain this little (but systematic) difference: 

• Experimentally, wicking is reported as soon as the first grain detaches. Thus 
if the actual critical contact angle 6^ locally differs from our experimental 
wicking criterion will determine only the highest possible critical contact 
angle. 

• In Fig. |2]:, one can see that the diameter of the grains is not perfectly fixed, 
which corresponds to the standard deviations reported in table [T} A small 
polydispersity can locally modify the geometry inside the pile, and thus 
the local value of the critical contact angle 6^ can be different from the 
monodisperse one 6q. In section 4, an elementary case of polydispersity, 
two-sized piles, will be investigated. 

• The interface was supposed to be flat, which is not the case if the difference of 
pressure through the interface is not negligible. In section 5, we will discuss 
this effect of pressure in the case where the thickness of the pile becomes 
comparable to the capillary length. 

• The determination of 6q assumes a close packing of spherical particles. As 
discussed in section 6, if the compacity of the pile is lower, or if the particle 
are slightly elongated, defects in the pile can appear, which modifies the 
wetting transition. 

4 Polydispersity 

As observed in Fig. |2(c)[ there is a small dispersion of bead radius. Ban et al. 
considered a polydisperse packing and anticipated theoretically that impregnation 
should be easier than in a monodisperse pile, although they found no difference 
experimentally Il28ll . 

The theoretical determination of the critical angle is similar to the monodisperse 
case in 3 dimensions, but geometry is modified by the bidispersity of grains (Fig. 
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Figure 9: Limit case for the impregnation of a bidisperse pile of grains. The dashed 
circle represents the spheres of the lower layer out of the projection plane. 



|9]). Eq. ^ can now be written for the lower layer: 

S = ]^icos0. (9) 

In addition, we have: 

R^ + S= {Ri + R2) cos oc, (10) 

where oc is determined by: 

sma = =— . (11) 

Using equations (L^ and (H) , the value of the critical angle 6^, function of 
the ratio r = ^ < 1, is: 



l + 2r-% 



cos r = ^ . (12) 

r 

When the upper sphere is much larger than the lower ones (r 0), the critical 
contact angle goes to 0"". Interestingly, this equation remains the same if the upper 
bead is the smaller one (r > 1), and it is valid as long as the upper sphere is large 
enough to stand upon the lower layer, that is r < 3 + 2\/3. In the monodisperse 
limit (r ^ 1), it gives the same result than Eq. ([8]). Additionally, for a small 
polydispersity, this equation yields at leading order in (r — 1): 

cosr = cos0^ + (1 - ^)(r - 1). (13) 

Since 1 — ^ is negative, the critical contact angle increases with r and is higher 
than for a smaller sphere on a bigger one (r > 1). This situation is the one 
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(a) (h) Vf = 0.80, = 48° =b 5° (c) Vf = 0.97, = 36° =b 5° 



Figure 10: (a): Set-up of the bidisperse experiment, (b): The R2 = 256 ± 13 }im 
beads stay on top of the base monolayer of radius Ri = 52 ± 2 }im. (c): Same 
experiment, with a lower contact angle. Only the top poles of the R2 beads are 
visible. The R2 spheres pass partially through the monolayer. 



favorable to wicking and thus should determine the path followed during the 
wicking of a random bidisperse pile. That said, for applications where wicking 
has to be guaranteed or prevented, Eq. ([12]) allows identifying the most stringent 
contact angle condition, which has to be met in the most unfavorable case of a 
local segregation. 

In order to investigate the effect of polydispersity on impregnation, we used 
beads with two distinct sizes Ri and R2. Different configurations can be obtained, 
with various associated compacities. We will focus on two layers system, each 
layer formed of beads of a given size. 

A monolayer of beads of radius Ri is deposited on the surface of the bath, 
with a contact angle 6. Then beads of radius R2 are added one by one on the first 
layer (Fig. flOk). We observe two regimes, depending on 6.1f6> 6^, the beads of 
radius R2 stay upon the monolayer of beads of radius Ri, and no impregnation 
is observed (Fig. p^). Conversely, for 6 < 6^, the beads of radius R2 are wetted 
and pass partially through the interface, between the beads of the monolayer 
(Fig. p^). Some Ri beads detach from the interface and sink into the liquid while 
only the north pole of the upper sphere remains dry, which corresponds to the 
impregnation of a bidisperse pile. Combinations of spheres presented in table [l] 
are used to obtain different values of r between 0.2 and 2.5, and the corresponding 
critical contact angles 9^ are reported in Fig. 11 , for comparison to the model. The 
critical contact angle of wetting is indeed found to be a function of the ratio of 
radii r — ^ 



R 



2 



Eq. ([12]) is coherent with experimental data, showing the effect of polydispersity 
on the wicking of elementary piles is well captured by our model, even if the value 
of the critical contact angle seems slightly underestimated (Fig. [TT]). This effect 
can be due to small defects in the packing of the base monolayer, that will be 
investigated in section 6. However, this experiment emphasizes again how critical 
is the role of geometry in the wicking of grains. 

Coming back to the "monodisperse" experiments of section 1, the size of beads 
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Figure 11: Wicking phase diagram for a bidisperse powder, in terms of contact 
angle 6 versus the ratio of radii of the spheres ^ = Dots are experimental data 
of the critical contact angle 9^ against r and the solid line is Eq. ([12]). 



is slightly varies from a layer to another, hence modifying the local critical contact 
angle. The dispersion of radius typically yields a ratio r ^ 1.1 ([l]), which generates 
a critical contact angle 6^ of 52°: polydispersity can explain part of the difference 
between the experimental data in [l] and the theoretical value of 6q given by Eq. ([8]). 

5 Hydrostatic forcing of the wetting transition 




(a) /z < /z^ (b) /z = /z^ (c) /z > /z^ 



Figure 12: Forced impregnation experiments at an ethanol volume fraction Vf = 
0.33, varying the pile thickness h. The corresponding contact angle is 6 = 80° ± 
5"^ > . When h is bigger than h^, impregnation is observed. 
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Even if the thickness h of the dry pile was set as small as possible in section 
2, spheres can be coarse enough to make h comparable to the capillary length 

a = In this case, the hydrostatic pressure plays a role, changing the shape of 

the interface between beads and thus the limit for impregnation. In the previous 
section, we implicitly assumed /z = 0. The following experiment shows the 
influence of the pile thickness on impregnation. On a liquid with a contact angle 
6 > 6q (no wicking), we add beads to the surface to increase the thickness. As in 
section 2, we observe a wicking transition: for h lower than a value h^, the pile 
remains dry (Fig. 12 a). As the thickness h reaches h^, a few grains detach from the 
interface (Fig. [T2|3). For thicker piles, the dry powder is impregnated, and many 
beads fall in the bath (Fig. 12:). The reference /z = is taken at the lowest point of 
the meniscus in the transparent cell, before any grain is added. Fig. 13 shows a 
phase diagram in terms of contact angle and thickness, separating the "wet" state 
(wicking) from the "dry" one (non-impregnation). At a given height /z, wicking 
only occurs if 6 is smaller than a critical value 6^, which increases with h. The two 
domains are separated by a line, of slope 30 mm/rad for R = 52 }im . 

If the grain radius is comparable to a, or if there is a difference of pressure AP 
between the liquid and the gas, the interface gets curved between the beads, as 
illustrated in Fig. [14| The associated Laplace pressure compensates the pressure 
AP, and changes the condition for impregnation. The Laplace equation gives the 
radius of curvature of the meniscus at equilibrium where 7 is the surface 




Figure 13: Phase diagram in terms of the thickness h against contact angle Q. The 
dots correspond to the limit between wet and dry states and thus represents the 
critical thickness for contact angles 6 higher than 0q, for glass beads with radius 
R = 52 }im. The line is the linear regression with slope 30 mm/ rad . 
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(a) 



(b) 



Figure 14: Limit case for impregnation with forcing, and close up on the meniscus. 
Due to the curvature of the meniscus, the liquid reaches the second layer of grains 
while no impregnation would have been observed with a flat meniscus and the 
same contact angle. 



tension of the liquid. The sign of this curvature is related to the sign of the forcing: 
a larger pressure in the liquid helps the impregnation process, elevating the highest 
point of the meniscus (Fig. [14]). The expression of 5 changes, while Eq. ^ remains 
correct. If jS is the angle between the meniscus at the contact point with the sphere 



and the equatorial plane (Fig. 14), one can get: 



6 = Rcos{d - iS) + - cos |6). (14) 

The distance between the two contact points can be written as a function of 
either ^ or 6 — ^ (Fig. [14]), which leads to: 

^siniS = R(l-sin(e-^)). (15) 

To determine j6, we can consider a situation close to the flat meniscus. Assum- 
ing jS <C 1, and a radius of curvature of the meniscus larger than the diameter of 

7 



the grains (^^ <C 1), Eq. (L5^ simplifies at the leading order into: 



Hence Eq. ([T4| becomes: 



^ = ^(l-sin^). (16) 



/?2 A p 

S = Rcos0+— — cos^e. (17) 
2j 



As shown in section 3, the critical angle for impregnation with a flat meniscus is 
given by cos 9q = S/R. With a curved meniscus, the critical contact angle 6* is 
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0.125 




Figure 15: Normalized critical thickness of the powder as a function of the 
cosine of the critical contact angle cos 6^ for beads of radius R = 52 jim. 



the contact angle 6 which solves the second order equation ([T7|), which yields at 
leading order: 

cos e"" = cos - cos^ 0^. (18) 

If the pressure AP is negative, which means a higher pressure in the air, it is 
harder to impregnate the pile (Q^ < 6q): the liquid has to wet more the surface 
to wet the grains. Conversely, if the pressure is positive, a liquid with a contact 
angle higher than 9q can invade the pile. Note that in the case where the radius of 
the grains is comparable to the capillary length, the curvature is set by a balance 
between the weight of the sphere, buoyancy and surface tension. Its sign will 
depend on the balance of the first two forces: if the spheres are denser than the 
liquid, the curvature is positive, and piles of such spheres can be impregnated 
with a contact angle higher than 0q, as indeed observed for the larger beads in[lj 

As seen in Fig. [12| the pressure is set hydrostatically: we have AP = pgh, where 
h is the depth of the deepest point of the pile (Fig. p^). The cosine of the critical 
contact angle should thus be linear in h as deduced from Eq. ([18]). A liquid with 
a contact angle ^ > i^^vades the powder if h is higher than a critical value 
given by: 

2a^ cosg^-cosg 

R cos^e^, ' ^ ^ 

where cos^q = \J\ — 1 (Eq. ([8])) and ^ 3.1mm^ is almost constant in the 

range of ethanol volume fraction used for this experiment {a^ = 3.2 mm^ for 
Vf = 0.4, = 3.0 mm^ for Vf = 0.6). 
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In Fig. 15 the normalized critical height is observed to decrease with 



increasing cos 6 with a linear relationship, as predicted by Eq. ([19]). However, the 
experimental value of the slope is —0.6, an order of magnitude smaller than the 
one predicted by our model, ^^"2^ ^ This discrepancy on the slope can be 
due to the large errors on both and contact angle, measured in a narrow domain 
of cos 6. 



6 Eccentricity and defects in the pile 



(a) e = 105°, R = 51^m (b) = 105°, R = 100 ^nl 




(c) e = 45°, R = bl}im (d) = 45°, R = 100 ^m 



Figure 16: Compacity of a monolayer of radius R on liquid interfaces with a contact 
angle 6. Pictures are taken from above through a binocular x20. 

If defects are present in the packing of the pile, they can affect the impregnation. 
Indeed, the observation of a monolayer of beads (Fig. [16]) shows that packing is not 
exactly a close packing. Packing depends on the way the monolayer is prepared 
experimentally, and defects appear especially if the liquid does not impregnate 
a monodisperse pile {6 > 6^). In this section, the effect of a gap between spheres 
will be discussed. 

In order to evaluate the gap, we compare the surface fraction occupied by the 
sphere to a compact situation. In a 2-D close packing, the surface fraction occupied 
by spheres is = (Fig. 17(a)). If the spheres have a gap of leR between 

them, the compaction decreases to O = j^^^ (Fig- 17(b)), so the gap ratio can be 
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(a) (b) (c) 

Figure 17: Surface fraction occupied by spheres at the liquid air interface: in the 
triangle, the surface of the grains is 12. In the close packing case (a), the total 
surface is \f?>V? while it is \/3(l + if there is a gap leR between spheres (b). 
(c): As model defects in the first layer of the pile, we consider an average gap leR 
between the grains. 



expressed as: 

^-1- (20) 




With pictures of packed spheres (Fig. 16), one can measure the observed 
compacity O and thus deduce the mean value of e reported in|2} 

Experimentally, the monolayer has less defects when made upon a wicking 
liquid for a monodisperse pile: it is then possible to add a small excess of grains in 
order to increase the packing. Note that in the polydisperse case, for r > 1 we are 
in the opposite situation, thus we expect the base monolayer to be poorly packed. 

We will focus on the situation favoring the impregnation, in particular when 
the defect is present in the first layer of grains. The defect is set to be a gap 7.eR 



between the beads of the first layer, as presented in the 2-D case in Fig. 17 (c). 

For a monodisperse pile, equations (|3]) and (|4]) are still correct. In the 3-D case, 
since the base of the tetrahedron is loose as presented in Fig. p7|(b), oc is given by: 

sin^ = ^(1 + e). (21) 

\j 

These equations lead to a modified critical contact angle, function of e\ 



cosr=^^(l-e-^)-l. (22) 

The comparison between experimental values (Table [T]) and the theoretical 
value leads to estimations of e in the monodisperse case: we find e = 4% 
for R = 25 Jim and R = 52 jitn, e = 7% for R = 100 jim and e = 12% for 
R = 256 jim. These results are coherent with the fact that smaller spheres form 
more compact layers at the surface of liquid, due to a relatively stronger interaction 
of the meniscus (Fig. 16(c) and Fig. |16(d)[ ) |[33l. These estimations are higher than 



the values measured on monolayers at the surface (Table [2]), but the estimated 
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R [fim] 


6 = 45° < r 


6 = 105° > e* 


52 ±2 


<D = 0.90, e = 0.4% 


<D = 0.72, e = 12% 


100 ±5 


<D = 0.87, e = 2.1% 


<D = 0.76, e = 9.2% 


256 ± 13 


<D = 0.86, e = 2.7% 


<D = 0.82, e = 5.2% 



Table 2: Estimated values of the surface fraction of beads O and average relative 
gap €, corresponding to the different situations and radius. 



values of e on monolayer correspond to mean values, whereas the experimental 
protocol is sensitive to local larger values of e on the surface. 

With the same arguments and those exposed in section 4, one can write the 
critical contact angle for a bidisperse pile with a defect 2eRi in the first layer: 



^l + 2r-^-¥(e+f)-l 

cos 0^ = ^ — , (23) 

r 

where r = ^ is the ratio of radii between the first and second layers, as defined 
in section 4. In the compact case {e 0), this equation gives the same result as 
Eq. ([12]). In Fig. [TlJ a systematic shift from the theory is observed, which can 
be partially explained with Eq. (|23]) with e = 0.6%, a value coherent with the 
observations from Fig. [I6j In order to test further this result, a precise control of 
the packing of the monolayer at the interface is needed, which exceeded the scope 
of our work. 



Conclusion 

In order to study the wicking in granular media, we carried out experiments 
on piles of glass beads at the surface of a bath, and shown the existence of a 
critical wetting angle 6^ below which wicking spontaneously occurs. This angle is 
significantly smaller than the one observed in capillary tubes (90^), and its value is 
close to the value expected from models, around 51^. The critical contact angle 
can be modified by the existence of a gradient of pressure through the liquid-air 
interface, or by the polydispersity of the grains. These effects can be used in 
industrial processes in order to either help or prevent the wicking of a powder, 
depending on the field of application. Moreover, this study emphasizes the crucial 
role of the geometry of the pores in the wicking of porous media. Focusing on this 
aspect, it may be possible to design the geometry of a porous media that will be 
spontaneously invaded by any liquid in contact (textile fibers), or one for which 
wicking only occurs for total wetting, as observed with super-oleophobic materials 



18 



References 

[I] Green, W. H.; Ampt, G. A. Studies on Soil Physics, Journ. of Agricultural Sc. 
1911, 4, 1-24. 

[2] Doerr, S. H.; Shakesby, R. A.; Walsh, R. R D. Soil water repellency: its causes, 
characteristics and hydro-geomorphological significance. Earth-Sc. Rev. 2000, 
51, 33-65. 

[3] Doerr, S. H.; Shakesby, R. A.; Dekker, L. W.; Ritsema, C. J. Occurrence, pre- 
diction and hydrological effects of water repellency amongst major soil and 
land-use types in a humid temperate climate. Eur. Journ. of Soil Sc. 2006, 57, 
741-754. 

[4] Amro, M. M. Factors Affecting Chemical Remediation of Oil Contaminated 
Water-Wetted Soil. Chem. Engineering & Tech. 2004, 27 890-894. 

[5] Taiz L.; Zeiger E., Plant Physiology, Sinauer 2010. 

[6] Ponomarenko, A.; Quere, D.; Clanet, C. A universal law for capillary rise in 
corners. Journ. Fluid Mech. 2011, 666, 146-154. 

[7] Hall, C; Hoff, W. D. Water Transport in Brick, Stone and Concrete, Second 
Edition, Taylor & Francis, London and New York 2012. 

[8] Staples, T. L.; Shaffer, D. G. Wicking flow in irregular capillaries. Colloids and 
Surfaces A 2002, 204, 239-250. 

[9] Brielles, N.; Chantraine, R; Viana, M.; Chulia, D.; Branlard, P.; Rubinstenn, G.; 
Lequeux, R; Mondain-Monval, O. Dissolution of a surfactant-containing active 
porous material. Journal of Colloid and Interface Science 2008, 328, 344-352. 

[10] Bico J.; Quere, D. Precursors of impregnation, Europhys. Lett. 2003, 61 348. 

[II] Alava, M.; Niskanen, K. The physics of paper. Rep. Prog. Phys. 2006, 69, 669. 

[12] Wall, P. E. Thin layer chromatography: a modern practical approach. Royal 
Soc. of Chem. 2005. 

[13] Bouaidat, S.; Hansen, O.; Bruus, H.; Berendsen, C; Bau-Madsen, N. K.; Thom- 
sen, P. Surface-directed capillary system; theory, experiments and applications. 
Lab Chip 2005, 5, 827-836. 

[14] Forny, L.; Marabi, A.; Palzer, S. Wetting, disintegration and dissolution of 
agglomerated water soluble powders. Powder Technology 2011, 206, 72-78. 

[15] Concus P.; Finn, R. On the behavior of a capillary surface in a wedge Proc. 
Nat. Ac. Sc. USA 1969, 63, 292-299. 

[16] de Gennes, P.-G.; Brochard-Wyart, R; Quere, D. Capillarity and wetting 
phenomena: drops, bubbles, pearls, waves. Springer 2004. 



19 



[17] Blow, M. L.; Kusumaatmaja, H.; Yeomans, J. M. Imbibition through an array 
of triangular posts, /. Phys.: Condens. Matter 2009, 22, 464125. 

[18] Hauksbee, F. R. S. Several Experiments Touching the Seeming Spontaneous 
Ascent of Water, Phil Trans, of the Royal Soc. 1708/1709, 26, 258-266. 

[19] Jurin, J. An account of some experiments shown before the Royal Society; with 
an enquiry into the cause of the ascent and suspension of water in capillary. 
Phil Trans, of the Royal Soc. 1717/1719, 30, 739-747. 

[20] Laplace, R S. Theory of Capillary Attraction. Supplements to the 10th book of 
Celestial Mechanics 1806/1807. Reprinted by Chelsea, New York 1966. 

[21] Lucas, R. Ueber das Zeitgesetz des kapillaren Aufstiegs von Fliissigkeiten, 
Kolloid-Z. 1918, 23, 15-22. 

[22] Washburn, E. W. The Dynamics of Capillary Flow, Phys. Rev. 1921, 17, 273-283. 

[23] Fries, N.; Dreyer, M. An analytic solution of capillary rise restrained by 
gravity. Journal of Colloid and Interface Science 2008, 320, 259-263. 

[24] Reyssat, M.; Sangne, L. Y.; van Nierop, E. A.; Stone, H. A. Imbibition in 
layered systems of packed beads, EPL 2009, 86, 56002. 

[25] Tsori, Y. Discontinuous Liquid Rise in Capillaries with Varying Cross-Sections, 
Langmuir 2006, 22, 8860-8863. 

[26] Czachor, H. Modeling the effect of pore structure and wetting angles on 
capillary rise in soils having different wettabilities. Journal of Hydrology 2006, 
328, 604-613. 

[27] Lago, M.; Araujo, M. Capillary rise in porous media. Physica A : Stat. Mech. 
and its App., 2001, 289, 1-17. 

[28] Ban, S.; Wolfram, E.; Rohrsetzer, S. The Condition of Starting of Liquid 
Imbibition in Powders, Colloids and Surfaces 1987, 22, 301-309. 

[29] Shirtcliffe,N. J.; McHale, G.; Newton, M. I.; Pyatt, F. B.; Doerr, S. H. Critical 
conditions for the wetting of soils. Applied Physics Letters 2006, 89, 094101. 

[30] Qian, B.; Shen, Z. Fabrication of superhydrophobic surfaces by dislocation- 
selective chemical etching on aluminum, copper, and zinc substrates. Langmuir 
2005, 21, 9007- 9009. 

[31] Specific gravity of mixtures of ethyl alcohol and water by volume and by 
weight. Handbook of Chemistry and Physics, Thirty-first edition 1949 1693-1699. 

[32] Keller, J. B. Surface tension force on a partly submerged body. Phys. Fluids 
1998, 10, 3009-3010. 



20 



[33] Mansfield, E. H.; Sepangi, H. R.; Eastwood, E. A. Equilibrium and mutual 
attraction or repulsion of objects supported by surface tension, Phil. Trans. R. 
Soc. Lond. A 1997, 355, 869-919. 

[34] Tuteja, A.; Choi, W.; Ma, M.; Mabry, J. M.; Mazzella, S. A.; Rutledge, G. C; 
McKinley, G. H.; Cohen, R. E. Designing Superoleophobic Surfaces, Science 
2007, 318, 1618-1622. 



21 



